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ABSTRACT 

Context. The origins of SO galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on 
environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a 
method for distinguishing between them. 

Aims. We aim to study a sample of six SO galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their 
stellar populations out to very large radii, to determine their kinematics in order to understand their origins. 

Mettiods. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these 
six systems * . 

Results. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what 
would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those 
derived at smaller radii from conventional spectroscopy. 

Conclusions. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of 
SO galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We 
show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most 
cases. 

Key words, galaxies: elliptical and lenticular - galaxies: evolution - galaxies: formation - galaxies: individual (NGC 1023, 
NGC 2768, NGC 3115, NGC 3384, NGC 3489, NGC 7457) - galaxies: kinematics and dynamics - galaxies: stellar content 



1. Introduction 

SO or lenticular galaxies seem to occupy a rather important 
place in the pantheon of galaxy types. Their smooth appear- 
ance coupled with their disk morphologies led iHubblel (119361) 
to place them at the centre of his sequence of galaxies between 



* (Tables 3 7 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsw eb.u-strasbg.fr/cgi-bin/qc at7J/A-l-A/) 



the smooth ellipticals and the disk-dominated spirals, suggesting 
that they may provide a key clue to understan ding the relation- 
ship between these galaxies. Further, 'Dressler'('l980) established 
that in all dense environments SO galaxies are the most numerous 
luminous systems, outnumbering even the ellipticals in the cen- 
tres of clusters, so clearly any theory of galaxy evolution has to 
contain a convincing explanation for the origins of SO galaxies. 

It is therefore rather frustrating that there is no consensus 
as to how these systems form. The lack of recent star forma- 
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tions implies that they have not had ready access to a sup- 
ply of cold gas in their recent history, but the reason for this 
lack of raw material remains a matter of debate. The quench- 
ing of star formation coul d be the result of galaxy minor merg- 
ers (iBoumaud et al.ll2005l) . analogous to the manner that more 
major mergers are believed to shut down star formation in ellip- 
tical galaxies. Alternatively, a normal spiral galaxy could sim- 
ply have had its gas stripped out by any of a range of more 
gentl e processes such as ram-pressure stripping or tidal interac- 
tions taunn & Gotl!'l972':'Arag6n-Salamanca'2008^. iOuilis et alJ 
[2000; Kronbergeret al. 2008; Byrd & Valtonen 199Cj). Thus, the 
most basic question of whether SO galaxies are more closely re- 
lated to ellipticals or spirals remains unanswered, as does the 
issue of whether multiple routes might lead to their formation 
depending on their surroundings. 

Fortunately, there are as-yet largely untapped clues to SO for- 
mation contained in their stellar kinematics. The simpler gas 
stripping processes would not have much impact on the colli- 
sionless stellar component, so one would expect an SO formed 
in this way to have the same stellar kinematics as its progen- 
itor spiral (Aragon-Salamanca 2008; Williams et al. 2010). On 
the other hand, a minor merger will heat the kinematics of the 
stellar component as well, leading to disks wi th substantially 
more random motion than is found in spirals (iBournaud et al.l 
l2005 f). More generally, elliptical and spiral galaxies obey quite 
tight scaling relations that relate their kinematics to their pho- 
tometric properties, so by investigating where SOs fall on these 
relations we can obtain a simple model-independent insight as 
to how closely they are related to the other galaxy types, and 
whether those connections vary with environment. 

The main problem in exploiting this information is that it is 
observationally quite challenging to extract. The outer regions 
of SO disks are very faint, so conventional absorption-line spec- 
troscopy struggles to obtain data of sufficient depth to derive 
reliable disk kinematics. Further, the composite nature of these 
systems, with varying contributions to the kinematics from disk 
and bulge at all radii, means that very high quality spectral data 
would be required to separate out the components reliably. 

Fortunately, planetary nebulae (PNe) offer a simple alter- 
native tracer of stellar kinematics out to very large radii. The 
strong emission from the [O iii] 5007 A line in PNe makes them 
easy to identify even at very low surface brightnesses, and 
also to measure their line-of-sight velocities from the Doppler 
shift in the line. To carry out both identification and kine- 
matic measurement in a single efficient step, we designed a 
special -purpose instrument, t he Planetary Nebula Spectrograph 
[PN.S; iDouglas et alJ (l2002h l. as a visitor instrument for the 
4.2m William Herschel Telescope (WHT) of the Isaac Newton 
Group in La Palma. Although o riginally i ntended to study el- 
liptical galaxi es (Romanowskv e t alJl2003Upouglas et alJl2007t 
ICoccato et al'.li2009MNapohtano et al.ll20Ml201 Ih . the"PN.S has 
already prov ed to be very effectiv e in exploring the disks of spi- 
ral galaxies (iMerrett et al.ll2006l) . so extending its work to SOs 
was an obvious next step. 

As a pilot study to test this possibility, we obs erved the 
archetypal SO galaxy NGC 1023 (Cortes i et aP 1201 iV and as- 
certained that it is possible to separate the more complex com- 
posite kinematics in these systems using a maximum likelihood 
method, which was sh own to give more r eliable results than a 
conventional analysis jNoordermeer et aLl 2008). Based on the 
success of this proof-of-concept, we decided to carry out a more 
extensive PN.S study of 6 SO galaxies that reside in a variety 
of environments, to see if this factor drives any variations be- 
tween them, and have obtained observations with PN.S of each 



galaxy of sufficient depth to obtain samples of typically 100 PNe 
in each. 

In this paper we present the data from the survey and assess 
its reliability as a kinematic tracer in all such systems. Indeed, 
one concern raised with using PNe to measure stellar kinemat- 
ics is that our understanding of the progenitors of these ob- 
jects is somewhat incomplete, so it is possible that they do not 
form a representative sample of the underlying stellar population 
(Dekel et al. 2005). This concern has been largely addressed for 
elliptical systems, where the PNe follow the same density pro- 
file and kinematics as the old stellar population where they over- 
lap, and have consistent specific frequencies implyi ng that they 
are not some random scattering of younger stars (Coccat o et al.l 
2009). However, SOs are a different class of system, so this is- 
sue needs to be looked at for these objects as well, to check that 
the concern can be similarly dismissed. The PN data sets are 
also more broadly useful to the galaxy formation commun ity for 
studying the outer regions of SOs (e.g. lForbes etani2012l) . 

Accordingly, this paper presents a description of the sample 
itself (Section |2|i and the resulting catalogues (Section |3]l. We 
also carry out the necessary checks on the viability of PNe as 
tracers of the stellar population by comparing their spatial distri- 
bution to that of the starlight, and, quantifying the normalisation 
of this comparison, the specific frequency of PNe (Section |4}. 
We present the resulting basic kinematics of the galaxies here 
(Section|5]) and draw some initial conclusions (Section|6]l. In the 
interests of clarity and the timeliness of the data, we defer the 
detailed modelling of the data for the full sample, together with 
its interpretation in terms of our understanding of SO galaxy for- 
mation, to a subsequent paper. 



2. The SO sample 

This project aims to study a sample of SO galaxies that spans 
a range of environments, to investigate whether their proper- 
ties vary systematically with their surroundings, as would be ex- 
pected if there are multiple environmentally-driven channels for 
SO formation. As usual in astronomy, we are faced by various 
conflicting pressures in selecting such a sample. First, in order 
to obtain reasonable numbers of PNe, the galaxies have to be 
fairly nearby, but this means that we cannot probe the rare dens- 
est environments of rich clusters. Accordingly, we compromise 
by studying a sample that spans the widest range available to us, 
from isolated SOs to poor groups to rich groups. Second, the req- 
uisite observations take a long time even with an instrument as 
efficient as PN.S, so the sample has to be quite modest in size, 
but there is also likely to be intrinsic variance due to the personal 
life histories of particular galaxies, which argues for as large a 
sample as possible. Accordingly, we select two galaxies from 
each environment to at least begin to explore such variations, 
giving a final sample size of six galaxies. 

Driven by these considerations, together with the observa- 
tional constraints of the William Herschel Telescope, the objects 
observed for this project were NGC 3115 and NGC 7457 as 
isolated SOs, NGC 1023 and NGC 2768 in poor groups, and 
NGC 3384 and NGC 3489 which are both in the richer Leo 
Group. The basic properties of these galaxies are summarised in 
Table[Tl while the details of the PN.S observations of the objects 
are provided in Table |2] 
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Table 1. Basic observational properties of the galaxy sample. 



Name 


Type 


D 


cz 


Ngroup 


K 






[Mpc] 


[km/s] 




[mag] 


NGC3115 


SO-edge-on 


9.68 


663 


1 


5.88 ±0.017 


NGC 7457 


SAO-(rs) 


13.2 


812 


2 


8.192 ±0.024 


NGC 2768 


E6 


22.4 


1373 


3 


6.997 ±0.031 


NGC 1023 


SBO-(rs) 


11.4 


637 


8 


6.238 ±0.021 


NGC 3489 


SABO+(rs) 


12.1 


677 


14 


7.37 ± 0.009 


NGC 3384 


SBO-(s) 


11.6 


704 


14 


6.75 ±0.019 



Notes: Col. 1: Galaxy name. Col. 2: Morphological classification (NED). Col. 3: Distance to the galaxy iTonrv et al.ll2001h. C ol. 4: Systemic 
velocity (NED). Col. 5: Number of members in the group the galaxy belongs to iFouaue et al.ll992l : lde Vaucouleurs et al.ll976l) . Col. 6: K-band 
apparent magnitudes (2MASS). 

Table 2. Basic observational properties of the galaxy sample. 



Name 


date 


Texp 


seeing 


filter 


A 


Vmin 


Vc 


Vmax 


Npn 






[hrs] 


[arcsec] 




[A] 


[km/s] 


[km/s] 


[km/s] 




NGC 3115 


2007 - 2010 


10.6 


1.2 


ABO 


5026.00 


-168 


1150 


2467 


189 


NGC 7457 


2001 


12.05 


1.15 


B6 


5026.74 


258 


1194 


2129 


113 


NGC 2768 


2007 


15.0 


1.4 


BO 


5033.90 


686 


1623 


2560 


315 


NGC 1023 


2007 


3/2.7 


1.0 


ABO 


5026.00 


-168 


1150 


2467 


204 


NGC 3489 


2007 - 2011 


7.8 


1.4 


ABO 


5026.00 


-168 


1150 


2467 


60 


NGC 3384 


2003 -2011 


5.95 


1.6 


ABO 


5026.00 


-168 


1150 


2467 


93 



Notes: Col. 1: Galaxy name. Col. 2: Years in which the observations were carried out. Col. 3: Effective exposure time (exposure time normalised 
at 1" seeing). NGC 1023 has been observed in two fields of view and the tabled exposure times correspond to the East and West fields 
respectively. Col. 4: Average seeing. Col. 5: Filter used. Col. 6: central wavelength of the filter. Col. 7/9: corresponding velocity window. Col. 
8: central velocity. Col. 10: number of observed PNe. 



2.1. Companion galaxies 

Since some possible channels for transformation into SOs in- 
volve interactions with companion galaxies, we now examine the 
cases where sample galaxies have relatively close neighbours. 

NGC 1023 has a low surface-brightness compan- 
ion, NGC 1023 A, w hich lies close to the main system 
('Capaccioli et al."l986). In their catalogue of PNe in NGC 2013, 
iNoordermeeret al. (2008) found that 20 of the 204 detected 
PNe were in fact associated with NGC 20 13 A. 

NGC 3115 is associated with MGC-1-26-12 (NGC 3115B) 
and MGC- 1-26-21, and it is surrounded by some Hi clouds. 
MGC-1-26-12 and MGC- 1-26-21 are two faint galaxies with I- 
band and /T-band apparent magnitudes of 14.26 mag and 10.05 
mag respectively (Doyle et al. 2005). MGC 1-26-21 is shown in 
the left panel of Figure[T] its distance from the main galaxy is 17' 
and it has a diameter of 70.3" Q MGC-1-26-12 is too far away 
from the central galaxy to be seen in Figure JI] being at a distance 
of 43 '. Its major axis measures 12.10" (Ide Vaucouleurs et al.l 
Il991h . 

NGC 7457 i s acc ompanied by UGC 12311 
(Ide Vaucouleurs et al.l Il976 h. which is an irregular galaxy 
with a A'-band apparent magnitude of 12.4, and a diameter 
of 23.6", and it lies at a distance of 5.7' from the centre of 
NGC 7457. It is shown in the right panel of Figure[T] 

With the possible exception of NGC 1023 A, all of these 
companions are faint and relatively distant, so are unlikely to 
have played a significant role in the evolution of the sample 
galaxies. 



' (2MASS , ' http://irsa.ipac.caltech.edu/) 



3. Data reduction and catalogues 

The data from PN.S are rather unusual in form, and require a 
custom-built analysis package for their analysis. The details of 
the pipeline are pre sented in Douglas et al. (2002, 2007) and 
Me rrett et al.l (l2006h : here we give a brief summary of the main 
steps, describe some recent upgrades to its capabilities, and tab- 
ulate the resulting catalogues of PNe. 

3. 1 . The Planetary Nebula Spectrograph and the reduction 
pipeline 

The PN.S, described in detail in lDouglas et al.l (l2002 l). is a slit- 
less spectrogra ph that works using a cou nter-dispersed imagin- 
ing technique (iFehrenbachI 1 1 948 . 1947a'b). Two images of the 
same field are taken at the same time through a slitless spec- 
trograph equipped with a narrow -band [OIII] filter to restrict 
the spectral range, and they are dispersed in opposite directions. 
Thus, stars in the field appear in the image as short segments of 
spectrum around 5007 A while the PNe appear as point sources 
due to their monochromatic emission. The exact position of a 
PN [Oiii] 5007 A line in each image depends on both its po- 
sition on the sky and its observed wavelength, so matching up 
the detected PNe in pairs from the images dispersed in opposite 
directions makes it possible to obtain both their positions and 
velocities in a single step. 

To extract this information, we have developed a specialised 
reduction pipeline for the instrument. The first step in the reduc- 
tion is to create individual distortion-corrected images. This step 
includes the creation of a bad pixel mask, cosmic ray removal, 
and finally the undistortion of each image using arc lamp expo- 
sures obtained through a calibrated giid mask. As a second step, 
all images for a given galaxy are aligned and stacked, weighted 
according to the seeing, transparency etc. Although stars appear 
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Fig. 1. ESO Digital Sky Survey (DSS) images of NGC 3115 {left panel, 60' x 60') and NGC 7457 {right panel 14.8' x 14.8') and 
their companion galaxies, which are marked with red circles. North is up, and East to the left in both panels. 



elongated into trails in the dispersed images, their centroids are 
still sufficiently well defined for them to be used to perform an 
astrometric calibration, based on the corresponding star coordi- 
nates in the USNO catalogue. At this point, the images are fully 
reduced, and we are ready to extract the PNe data from them. 

The procedure for extracting a catalogue of PNe from 
these images is semi-automated, initially using Sextractor 
dBertin & Arnouts 1996) to find potential objects. As explained 
above, PNe in the images appear as point-like sources, since they 
emit only in the 5007 A line within the wavelength range se- 
lected by the filter. With Sextractor, experimentation shows that 
these real sources can be distinguished from residual cosmic rays 
and other image defects by selecting only those objects that com- 
prise a group of at least 5 contiguous pixels whose individual 
fluxes are at least 0.7 times the rms of the background. 

There is, however, still the possibility of contamination from 
foreground stars and background galaxies. Fortunately, such 
sources are elongated along the dispersion direction due to their 
continuum light emission, with the length of elongation dictated 
by the filter bandpass. To identify such objects automatically, we 
create two further images. The first is an unsharp mask formed 
by smoothing the original spectral image with a 49 x 49 median 
filter, and subtracting the smoothed image from the original. The 
second is obtained by convolving this median-subtracted image 
with an elliptical Gaussian function elongated in the dispersion 
direction. The resulting images are almost the same for extended 
stellar trails, but the point-like PNe are fainter in the smoothed 
version. Thus, for either arm, we can extract a catalogue that 
excludes any continuum sources by matching up source detec- 
tions, and selecting only sources whose difference in magnitude 
between the median-subtracted image and its smoothed counter- 
part is larger than 0.5 mag (again tuned by experimentation). 

Although this automated process, which we have used in our 
previous pipeline processing, detects the vast majority of PNe, 
there are still a number of missed detections when the object is 
projected near a star trail. We have therefore developed a fur- 
ther technique to subtract out the star trails as an aid to detecting 
PNe, illustrated in Figure |2] We create a 25 x 1 median-filtered 



image (top left panel) and subtract it from the science image, 
obtaining the median-subtracted image shown in the top mid- 
dle panel. This image smooths away most of the star trails and 
leaves PNe almost unaffected. However, there are still residual 
features from the brighter parts of the brightest star trails that 
could be misidentified as PNe. We must therefore create a mask 
to exclude these regions. To do so, we further apply a 81 x 81 
median filter to the 25 x 1 median filtered image to obtain an 
image that picks out the stars while removing the PNe. To turn 
this image into a mask, we run Sextractor on it and perform a 
morphological dilation on the segmentation map shown in the 
bottom left panel, by blurring (bottom middle panel) and thresh- 
olding (bottom right panel) the image. Applying this mask, by 
replacing masked values with zero, then yields the final trail- 
free data frame shown in the top right panel, which can then be 
pipeline processed in the way described above, with no residual 
confusion due to foreground stars. 

With this modification, the only PNe missed by the auto- 
mated process are those that lie very near the edge of the image, 
where the effectiveness of median filtering is limited, or near 
the centre of the target galaxy, where the bright continuum is 
noisy and difficult to subtract fully. However, the numbers are 
now sufficiently small that any remaining PNe are easily added 
in manually. 

Finally the catalogues obtained in this way for the two arms 
of the spectrograph are matched up. Any pairs that appear at the 
same location in the undispersed direction to within two pix- 
els, so have consistent spatial coordinate, and lie within a few 
hundred pixels of each other in the dispersed direction, so are 
consistent with the likely offset due to their Doppler shifts in 
the counter-dispersed frames, are identified as images of a sin- 
gle PN. Any object that does not appear in both arms is rejected 
as spurious. As a final check, each putative detection is assessed 
manually by three independent team members, by checking its 
location in the original stacked images, and an object is accepted 
in the catalogue if at least two classifiers confirm that it is real. 
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Fig. 2. Different steps of the star trail removal routine, the top right panel show the resulting image. Refer to text for major details. 



The wavelength and astrometric solutions derived above then 
allow us to obtain the mapping 

{xi,yi,Xr,yA ^ {xQ,yo,A}, (1) 

where xi,yi,Xr,yr are the coordinates of the PNe in the left and 
right-dispersed images, xq and yo are the true sky coordinates 
and A is the wavelength of emission for each object. The emis- 
sion wavelength translates directly into a line-of-sight velocity 
via the Doppler formula. 

EiTors are assigned to these velocities by propagating the er- 
rors in the assigned values of xi,yi,Xr,yr through the transfor- 
mation, further allowing for the uncertainty in the dispersion so- 
lution itself. A check on these uncertainties is provided by the 



fact that the observation of NGC 3384 presented here overlaps 
with th e field containing NGC 3379 published by Douglas et alj 
Matching the catalogues results in 25 objects in both 
observations for which the RMS difference between the data is 
20 km s"' , implying an average eiTor in each individual observa- 
tion of 14 km s ' , consistent with what the internal error analysis 
assigned. 



3.2. The catalogues 

The resulting catalogues of PNe are presented in Tables |3]-|7] The 
data for NGC 1023 are not shown here, since they were pub- 
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lished in iNoordermeer et al.l (l2008l) . For completeness, Table [T] 
also contains the 9 PNe identified as belonging to NGC 3384 in 
the adjacent observation of NGC 3379 by Douglas et al. (2007). 
As well as the positions and line-of-sight velocities of the PNe, 
these tables also show the instiTimental magnitude of each ob- 
ject, as derived from the images; since we only use these mag- 
nitudes to assess the strength of the detection and to investi- 
gate the shape of the planetary nebula luminosity function (see 
SectionUJi, no photometric calibration is required. 



4. The PNe Spatial Distribution and its 
Normalisation 

As discussed in the Introduction, one concern with using PNe as 
tracers is that they may not be representative of the whole stel- 
lar population, but rather a subcomponent, which may introduce 
a bias in the dynamical analysis. This issue has already been 
thoroughly addressed in the context of elUpticals (ICoccato et alJ 
120091; iNapolitano et al.ll2009l 1201 lb . but now we need to estab- 
lish the same point for SOs. Indeed, similar findings for SOs and 
ellipticals would not only confirm the effectiveness of the tech- 
nique, but might also allow us to draw connections between the 
two types of galaxy. 



4.1. Spatial distribution 

The simplest test we can make is to compare the spatial distribu- 
tion of PNe to that of the underlying stellar population: if PNe 
are effectively a random sampling of stars, then, to within a nor- 
malisation factor, the two should be the same. 

We have therefore fitted elliptical isophotes to images of the 
galaxies in the sample, and plotted the surface brightness profiles 
extracted from these fits in Figure |3] In most instances, the im- 
ages used were the A"-b and observations from the 2MASS sui-vey 
(ISkrutskie et al.ll2006l) . In the case of NGC 3489, the galaxy is 
too faint to be reliably fitted using this data, b ut fortunately thi s 
galaxy lies in the region surveyed by SDSS ( York et al.llBooh . 
so we used their z-band image as a substitute: since we are 
only interested in matching the shape of the brightness profile, 
this red band traces the old stellar population as well as the K- 
band, and no colour coiTection is required. Finally, in the case 
of NGC 1023 we made use of a deep /?-band image that we had 
previously analysed to determine the size of the region compro- 
mised by light from NGC 1023 A, to exclude it from the analysis 
(ICortesietal.ll201lh . 



Table 3. Catalogue of PNe in NGC 3115. 



Table 4. Catalogue of PNe in NGC 7457. 



ID 


RA 


Dec 


Vhelio 


m 




J2000 


J2000 


[km/s] 


[mag] 


PNS-EPN-N7457-001 


23:00:42.75 


30:10:10.7 


993.94 


20.4 


PNS-EPN-N7457-002 


23:00:48.31 


30:08:58.2 


870.64 


20.24 


PNS-EPN-N7457-003 


23:00:51.14 


30:09:19 


886.94 


19.5 


PNS-EPN-N7457-004 


23:00:51.42 


30:09:44.1 


939.94 


19.69 


PNS-EPN-N7457-112 


23:01:8.34 


30:08:24.6 


755.84 


19.98 


PNS-EPN-N7457-113 


23:01:10.11 


30:06:58.2 


724.74 


19.83 



Notes: Details on the velocity and positional errors are given in 
Sectior im Velocities, measured from the [O in] 5007 A line, are 
corrected to heliocentric. The full catalogue is available at the CDS. 



Table 5. Catalogue of PNe in NGC 2768. 



ID 


RA 
J2000 


Dec 
J2000 


Vhelio 
[km/s] 


m 
[mag] 


PNS-EPN-N2768-01 


09:10:56.13 


60:03:21.9 


1181.3 


20.86 


PNS-EPN-N2768-02 


09:10:59.59 


60:02:13.7 


1133.9 


21.86 


PNS-EPN-N2768-03 


09:10:59.86 


60:03:38.5 


1134.4 


21.05 


PNS-EPN-N2768-04 


09:11:1.34 


60:03:3.1 


1172.1 


20.49 


PNS-EPN-N2768-314 


09:12:15.08 


60:03:8.9 


1556.4 


20.36 


PNS-EPN-N2768-315 


09:12:15.22 


60:00:3.6 


1595.6 


21.08 



Notes: Details on the velocity and positional errors are given in 
Sectior im Velocities, measured from the [O m] 5007 A line, are 
corrected to heliocentric. The full catalogue is available at the CDS. 

Table 6. Catalogue of PNe in NGC 3489. 



ID 


RA 


Dec 


Vhelio 


m 




J2000 


12000 


[km/s] 


[mag] 


PNS-EPN-N3489-01 


11:00:09.1 


13:53:21.3 


660.52 


16.49 


PNS-EPN-N3489-02 


11:00:09.2 


13:55:26.3 


663.62 


16.59 


PNS-EPN-N3489-03 


11:00:09.6 


13:51:25.7 


1012.62 


19.17 


PNS-EPN-N3489-04 


11:00:10.4 


13:53:45.1 


1418.92 


16.49 


PNS-EPN-N3489-59 


11:00:23.7 


13:54:43.8 


742.12 


15.82 


PNS-EPN-N3489-60 


11:00:27.4 


13:54:41.6 


591.12 


16.62 



ID 


RA 


Dec 


Vhelio 


m 




12000 


J2000 


[km/s] 


[mag] 


PNS-EPN-N3489-01 


10:04:51 


-07:44:57 


645.58 




PNS-EPN-N3489-02 


10:04:54.9 


-07:43:55.5 


845.18 


18.96 


PNS-EPN-N3489-03 


10:04:55.5 


-07:43:37.4 


307.08 


19.1 


PNS-EPN-N3489-04 


10:04:57.1 


-07:46:22.8 


568.28 


18.65 


PNS-EPN-N3489-188 


10:05:36.2 


-07:41:28.1 


804.58 


17.56 


PNS-EPN-N3489-189 


10:05:36.5 


-07:41:49.7 


722.08 


19.17 



Notes: Details on the velocity and positional errors are given in 
Sectior im Velocities, measured from the [O m] 5007 A line, are 
corrected to heliocentric. The full catalogue is available at the CDS. 



Typically, the PNe extend to larger radii than the reliable 
■photometric fits shown in Figure |3] In order to extrapolate the 
comparison to these large radii, we have also fitted the images 
with a simple two -component disk-i-bulge model using GALFIT 



Notes: Details on the velocity and positional errors are given in 
Sectior [3Tn Velocities, measured from the [O iii] 5007 A line, are 
corrected to heliocentric. The full catalogue is available at the CDS. 



file, resulting in the lines shown in Figure |3] 

In order to compare these profiles to the PNe data, we need 
to correct the latter for incompleteness. Since PNe are primarily 
lost against the bright central parts of the galaxy, failing to cor- 
rect for this effect would grossly distort their spatial profile. We 
carry out t his correction i n two stages, following a process de- 
veloped bv lCoccato et al.l ( l2009h . First, we randomly place simu- 
lated PNe of varying magnitudes in model images that reproduce 



6 



Cortesi et al.: PN.S survey of SO galaxies 



Table 7. Catalogue of PNe in NGC 3384. 





ID 


RA 


Dec 


Vhelio 


m 






J2000 


J2000 


[km/s] 


[mag] 


PNS 


EPN-N3384-01 


10:47:54.8 


12:36:48.5 


822.6 


18.75 


PNS 


EPN-N3384-02 


10:47:54.9 


12:36:15.8 


1038.7 


19.01 


PNS 


EPN-N3384-03 


10:47:55.1 


12:36:16.8 


868.1 


18.98 


PNS 


EPN-N3384-04 


10:47:57.1 


12:36:02.8 


908.7 


18.86 


PNS 


EPN-N3384-94 


10:48:31.2 


12:39:02.6 


646.9 


18.41 


PNS 


EPN-N3384-95 


10:48:33.2 


12:36:26.9 


629 


18.23 



Notes: Details on the velocity and positional errors are given in 
Sectior fm Velocities, measured from the [O iii] 5007 A line, are 
corrected to heliocentric. The full catalogue is available at the CDS. 



s; 40 




Fig. 4. Recovered number of point like sources as a function of 
magnitude. The simulated PNe are added on an artificial image 
that mimic the sky and CCD noise, the recovered number de- 
creases for fainter objects, see text for further details. 



the noise properties of the real data, to establish the brightness 
level to which PNe can in general be reliably detected. We quan- 
tify this limits by mgo, the magnitude limit at which 80% of the 
simulated PNe are recovered. An illustration of this process is 
shown in Figure |4] Next, we populate the real science data with 
simulated PNe drawn from the known planetary nebula luminos- 
ity function (PNLF), brighter than an mgo completeness limit, 
and see what fraction are lost against the bright galaxy contin- 
uum, star trails, etc, as a function of position, as illustrated in 
Figure|5] Finally, we correct for this incompleteness in the num- 
ber counts of PNe as a function of radius, plotted in Figure |3] 

In the central parts of these galaxies, incompleteness is too 
great to be corrected reliably, but, as Figure |3] shows, in the re- 
gion of overlap with the photometric data, out into the extrap- 
olation of the model fit, the agreement between the profiles of 
PNe and the stellar continuum is very good, typically over at 
least a factor of ten in brightness. There is therefore no evidence 
that the PNe are anything other than a fair sample of the stellar 
population. 



4.2. The specific frequency of PNe 

Having established that the proportionality between PNe num- 
bers and stellar luminosity seems to be valid, we can now calcu- 



Fig. 5. Number of detected point like sources with m > mgo as a 
function of radius; the simulated objects are added to the science 
image. The incompleteness at small radii is caused by the bright 
continuum of the galaxy, see text for further details. 



Norm = 0.1518 

Cut-off =17.176t 
N(M) = 0J52e<'''='/'' [1 - a'"]'-''] 




insfrumenfal magnitudes 



Fig. 6. The observed PNLF and the fitted analytic function 
(ICiardulloetal.lll989h of NGC 3115, see text for further details. 



late the proportionality constant. This factor, usually designated 
a ( Jacob v 1980), is of interest because it relates to the stellar evo- 
lution of the particular population. So, for example, assessing 
how similar it is between ellipticals and SOs provides an inde- 
pendent route for determining how closely related these objects 
are, in this case in terms of the evolution of their stellar popu- 
lations. Such a similarity would also add further justification for 
the view that PNe in SOs have a similar relation to the underlying 
stellar population as those in ellipticals, and hence that they can 
be used as similar kinematic tracers. 

In fact, there is already evidence that this proportionality 
varies in interesting ways with the global properties of the 
galaxy. Although a is clos e to universal in value, it does vary 
systematically with colour ( Hui et alll993h . with the reddest el- 
lipticals being poorer in PNe per unit galaxy luminosity than spi- 
rals (Bu zzoni et al.i2006 ). This is somewhat counter-intuitive, as 
one might expect a star-forming spiral to have more of its lumi- 
nosity provided by m assive young stars that will never form PNe 
(iBuzzoni et alJl2006 ). However, there are other ways in which 
the a parameter might be reduced. For example, if PNe have 
systematically shorter lifetimes in elliptical galaxies, one would 
expect to see fewer of them in any snapshot view. Alternatively, 
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KGC 3115 



NGC 7457 




26 - 

aa : 



100 200 300 

R [arcsec] 

KGC 1023 





40 60 BO 

K [arcsec] 

NGC 2768 




KGC 3489 



NGC 3384 




60 80 100 

R [arcsec] 



100 150 
H [arcsec] 



Fig. 3. Surface density profiles of PNe and stars. Open diamonds: surface brightness profiles obtained with ellipse fitting of the 
galaxy image. Filled circles: PNe number density, computed as -l.SLogNcjA, where Nc is the completeness corrected number of 
PNe with m < m^o in the radial bin, and A is the area of the circular anulus. Blue lines represent the extrapolation of the sersic + 
exponential disk fit to the stellar surface brighness data. 



some stars with particular chemical properties might avoid be- 
coming PNe at all, perhaps by losing their envelopes at an ear- 
lier stage. Support for this suggestion comes from the rather tight 
anti-correlation observed between a and the amount of far ultra- 
violet (FUV) emission from the galaxy: in an old stellar popula- 
tion, FUV emission could well arise from AGB stars that have 
alrea dy lost their envelop es, thus avoiding later evolution into 
PNe (iBuzzoni et alJl2006l) . 

We therefore now calculate the a parameter for the sample 
SO galaxies, to see how they compare t o those already deter- 
mined for ellipticals (ICoccat o et al.' 2009). The value of a obvi- 
ously depends on how far down the PNLF one counts PNe away 
from the brightest magnitude observed. This brightest magnitude 
is well-defined by a sharp cut-off in the PNLF at a magnitude m*. 



which we can measure by fitting the lCiarduUo et al.l(ll989h ana- 
lytic approximation to the data, as illustrated in Figure |6](kept in 
instrumental magnitudes since an absolute value is not required). 
The value of a then depends on how far below m* one counts 
PNe. For a direct comparison with the e llipticals, we follow the 
same definition as ICoccato et all (l2009l ). which involves count- 
ing the number of PNe with apparent magnitudes down to m* -i- 1 . 
We then calculate aBj.o, which is the ratio of this number of PNe 
to the total B-band luminosity over the range of radii that PNe 
have been observed. The radial range used for each galaxy in 
this survey is listed in Table |8] 

In practice, our data generally probe to magnitudes signifi- 
cantly fainter than m* + I: Table [8] lists the values of m* - mgo 
which quantifies how many magnitudes below m* the data re- 
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main reasonably complete. In order to use all this information 
in determining a, we calculate the constant of proportional- 
ity counting PNe all the way down to this limiting magnitude, 
aB,m'~mm' ™d then use the known shape of the PNLF to correct 
its value to the adopted common metric quantity. 



dBA.O - 0!B«i'-n 



J2" F{m*,m)dm 

) r»i8o r,, I 77~, T' 
J,„. r(m ,m)am 



(2) 



where F(m*,m) is the universa l analy tic approximation to the 
PNLF found bv'C iardullo et alJ (1 19891) . The values for agj o de- 
termined in this way are listed in Table |8] 

Having calculated these a values, we can now compare 
th em to the valu e s dete rmined in the same way for ellipticals 
bv lCoccato eTaP (l2009h . Figure |7] shows the resulting compar- 
isons when combining a values with optical colours, FUV ex- 
cess and absolute magnitude. Plotted against optical colour, the 
SOs a values are indistinguishable from ellipticals, and seem to 
show the same trend described above that bluer galaxies have 
higher specific frequencies of PNe. Similarly, when we inves- 
tigate their FUV fluxes, obtained from the Galaxy Evolution 
Explorer (GALEXfl and corrected for extinction using the re- 
lation Afuv = ^316E{B - V) dWvder et al.ll2005h . we find that 
the SOs lie in the same part of the plot as ellipticals, although the 
range of FUV excesses for the SOs is too small to really see if 
they follow the same strong trend. Finally, we compare a to the 
B-band absolute magnitudes, which were obtained from the ap- 
parent magnitudes in the RC3 catalogue and the distance mod- 
uli of iTonrv et all ( 1200 lb . shifted by -0.16 magnitudes to take 
into a ccount the newer Cep heid zero-poin t of Freedm an et al.l 
(1200 1). for consistency with Coccat o et al.l (120091) . There is an 
interesting suggestion of a coiTelation between a and FUV - V 
for SOs in Figure |7] which is strengthened by the fact that two of 
the elliptical galaxies that lie closest to the SOs (NGC 4564 and 
NGC 3377) were found bv lCoccato eTa l. (2009) to have disk Uke 
kinematics, and thus may well be SOs themselves. However, the 
numbers are still sufficiently small that the only firm conclusion 
is that SOs lie in the same region of these plots as ellipticals, and 
hence that the stellar evolutionary properties that they probe are 
indistinguishable, adding further confidence to the use of PNe as 
a common steUar kinematic probe in these galaxy types. 



5. Velocity maps and kinematics 

In this section we present the resulting stellar kinematics. 
Figure |8] shows the PNe data superimposed on the Digitised 
Sky Survey (DSS) image of each galaxy in the sample. The in- 
completeness at small radii due to the bright galaxy continuum, 
discussed in Section |4l is very apparent. However, these plots 
also show the large radii to which PNe can be detected in these 
data, which underlines how nicely they complement conven- 
tional kinematic data: the point where the continuum becomes 
too faint to obtain good absorption-line data is just where PNe 
become readily detectable, and they remain present in significant 
numbers out to the very faint surface brightnesses illustrated in 
Figure |3] 

It is apparent from Figure |8] that most of these galaxies are 
characterised by obvious rotation, although there is also a sig- 
nificant amount of random motion. It also appears that the bal- 
ance between ordered and random motion differs from galaxy to 
galaxy; for example NGC 7457 seems to show more coherent 



GALEX data are available at: |http://galex.stsci.edu/GR6/| 



motions than NGC 3115. To quantify such differences in a pre- 
liminary kinematic analysis, we derive the rotation velocity and 
random motions assuming that each gala xy can be described as 
a single-component thin disk. Following ICortesi et al.l (l201 ih . a 
maximum likelihood analysis is employed to determine in a se- 
ries of radial bins the mean streaming velocity and the two com- 
ponents of random motion in the tangential and radial directions, 
here assumed to be coupled together by the epicycle approxima- 
tion. To minimise distortion of the fit by outlier PNe, the analysis 
iteratively rejects objects whose deviation from the mean veloc- 
ity in each bin is higher than 2.1 cr; in no case are more than 
7% of the total number of PNe rejected by this process, offering 
some reassurance that such a disk model is not unreasonable. 

As Figure |9]illustrates, this analysis confirms the impression 
that the balance between ordered and random motion varies sig- 
nificantly between galaxies. It also offers one final very useful 
check on the reliability of PNe as stellar kinematic tracers, as we 
can compare the results to those from conventional absorption- 
line spectra. Figure |9] therefore also shows such absorption -line 
kinematics for the major axes of these galaxies, as obtained from 
various publ ished sources ( Debattista et al. 2002; Caon et^ 
2000; Simie n & PmgnieJl997HNorris et al.ll2006tlFisherlll997ir 
Although the comparison is not exact, these systems are suffi- 
ciently close to edge on that the conventional major-axis velocity 
dispersion data should lie close to cr^ in the PNe data if both are 
measuring the same kinematic population. These plots underline 
the complementarity between conventional absorption-line data 
and PNe kinematics, as the latter take over just where the former 
run out. It is also clear that there is no inconsistency between 
these data sets, demonstrating once again that they are probing 
the same kinematic population. 

6. Conclusions 

This paper has presented the data for the first systematic sur- 
vey of SO kinematics as a function of environment, using PNe to 
trace their stellar velocities out to large radii. In the process, we 
have presented a series of tests on the viability of using PNe for 
such a study, demonstrating that they really are representative of 
the underlying stellar population. The PNe pass all these tests: 
they show the same spatial distribution as the starlight; their 
numbers are consistent with those found in comparable elliptical 
galaxies where their use as kinematic tracers is well established; 
and their kinematics join smoothly onto those derived from con- 
ventional stellar data at smaller radii. 

As well as demonstrating that PNe are good tracers of the 
stellar populations, and therefore that the catalogues presented 
in this paper are a unique dataset to study the kinematics and 
dynamics of lenticular galaxies, we have also shed some light 
on the nature of SO galaxies. In particular, the similarity of their 
a parameters to those found in elliptical galaxies implies that, 
at least at the stellar population level, these two classes of ob- 
ject are very similar. Unsurprisingly their kinematics are rather 
different, being in most cases dominated by rotational motion. 
There is however significant variety in these kinematics, with 
rotation velocities that vary in differing ways with radius and 
differing balances between ordered and random motion. The ini- 
tial inspection of these properties in Figure |9] does not reveal 
any obvious trends with the richness of environment, suggesting 
that this may not be a dominant factor in the choic e of chan- 
nel by which SOs form. However, Cortesi et al.l (1201 1 ) already 
established that these simple kinematic properties can be rather 
misleading due to cross-contamination between disk and bulge 
kinematics, and that a clear picture only emerges when the kine- 
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Table 8. ag^ix) parameter and UV excess for the galaxy sample. 



Name 


Npn 


NpNjCORR 


Rmin 


Rmax 


m, — mgo 


O^BA 


(FUV-V) 


Bt 








[arcs] 


[arcs] 


[mag] 


10' Lo 


[mag] 


[mag] 


NGC3115 


182 


186 


45.9 


389.8 


-3.0 


1.0 + 0.2 


7.6 + 0.1 


-20.03 + 0.13 


NGC 7457 


94 


156 


10.6 


115.4 


-1.9 


8.4 + 3.0 


8.1+0.1 


-18.69 + 0.32 


NGC 2768 


201 


289 


25.1 


328.0 


-1.6 


2.1+0.6 


7.6 + 0.1 


-20.89 + 0.34 


NGC 1023 


133 


236 


21.7 


332.4 


-2.5 


2.5+0.6 


7.6 + 0.1 


-19.10 + 0.28 


NGC 3489 


50 


75 


20.3 


146.0 


-3.3 


1.9 + 0.8 


6.6 + 0.1 


-20.05 + 0.22 


NGC 3384 


75 


77 


43.6 


228.1 


-2.6 


2.3 + 0.9 


7.5 + 0.1 


-19.41+0.19 



Notes: Col. 1: Galaxy name. Col. 2: Number of PNe with in < mso- Col. 3: Number of PNe with in < nigo once corrected for completeness. 
Col. 4: Galactocentric radius of the innermost detected PN. Col. 5: Galactocentric radius of the outermost detected PN. Col. 6: Difference 
between the PNLF bright cut-off (in,) and nigo (see Section|4j2j. Col. 7: a parameter in the i?-band, integrating the PNLF from in, to m, + 1 mag 
(see Section l42l l. Col. 8: UV excess, from the Galaxy Evol ution Explorer (GA LEX). Col. 9: S-band absolute magnitude, from the apparent 
magnitudes in the RC3 catalogue and the distance moduli of iTonrv et al.l ( 1200 ih . shifted by -0.16 magnitudes, see text for details. 

(B-V) [mag] 

0.4 0.6 0.8 1 1.2 

\ ^ \ 1 \ ^ 1 1 \ ^ 1 1 n ^ 1 \ \ ^ 1 

I E3 1 




5 6 7 8 9 -18 -19 -20 -21 -22 
(FUV-V) [mag] [mag] 

Fig. 7. Relationships for elliptical and SO galaxies between as^i.o parameter and UV excess (bottom left panel), absolute B-band 
m agnitude (bottom-righ t panel) and (B-V) colour (upper panel). Red open squares show the values of the elliptical sample presented 
in lCoccato et al.l ( l2009l) . while filled black circles correspond to SOs from the present study. 
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Fig. 8. DSS images of the sample galaxies, with the PNe positions marked with circles. The color of the circles represent the 
measured radial velocity of the PNe, according to the velocity scales below each panel. First row: field SOs, second row: SOs in 
groups, third row: SOs in rich groups. 
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Fig. 9. Kinematics of the SO sample derived using the disk model of lCortesi et"al] ( 1201 lb . First row: field galaxies, second row: 
group members, third row: rich group members. The black circles show the average velocity in the bin (whose size is represented by 
the horizontal error bars). The red squares and blue circles show the velocity dispersion along t he radial an d azimut hal directions re- 
spect i vely. The stellar ve l ocities and ve locity dispersions obtained from absorption line spectra (ICaon et al. 2000; Si mien & Prugnieil 
119971: iNorris et al.ll2006t lFishedll997i) are shown as empty black circles. The number of PNe left after the clipping performed in 
each UkeUhood fit is annotated at the top of the corresponding panel. 



matic components are carefully separated. Indeed, by breaking 
SO galaxies down into these components, we hope to find subtler 
indications as to where these systems lie in the overall story of 
galaxy morphology by, for example, comparing their spheroidal 
components to ellipticals and to the bulges of spirals. We will 
present such an analysis in a subsequent paper, to seek the defini- 
tive answers from the survey data presented here. 
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